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What is the form of the neutrino mass matrix which governs the oscillations of the atmospheric and 
solar neutrinos? Features of the data have led to a dominant viewpoint where the mass matrix has 
an ordered, regulated pattern, perhaps dictated by a flavor symmetry. We challenge this viewpoint, 
and demonstrate that the data are well accounted for by a neutrino mass matrix which appears to 
have random entries. 
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1 Neutrinos are the most poorly understood among 
known elementary particles, and have important con- 
sequences in particle and nuclear physics, astrophysics 
and cosmology. Special interests are devoted to neutrino 
oscillations, which, if they exist, imply physics beyond 
the standard model of particle physics, in particular 
neutrino masses. The SuperKamiokande data on the 
angular dependence of the atmospheric neutrino flux 
provides strong evidence for neutrino oscillations, with 
i/fj, disappearance via large, near maximal mixing, and 
Am 2 tm » 10~ 3 eV 2 (!]]. Several measurements of the 
solar neutrino flux can also be interpreted as neutrino 
oscillations, via v e disappearance g. While a variety 
of ArriQ and mixing angles fit the data, in most cases 
Am,Q is considerably lower than Am 2 tm , and even in the 
case of the large angle MSW solution, the data typically 
require AuIq « 0.lAm 2 tm The neutrino mass matrix 
apparently has an ordered, hierarchical form for the 
eigenvalues, even though it has a structure allowing large 
mixing angles. 

All attempts at explaining atmospheric and solar 
neutrino fluxes in terms of neutrino oscillations have 
resorted to some form of ordered, highly structured 
neutrino mass matrix Q. These structures take the 
form Mo + eMi + where the zeroth order mass 
matrix, Mo, contains the largest non-zero entries, but 
has many zero entries, while the first order correction 
terms, eMi, have their own definite texture, and are 
regulated in size by a small parameter e. Frequently 
the pattern of the zeroth order matrix is governed by 
a flavor symmetry, and the hierarchy of mass eigenvalues 
result from carefully-chosen, small, symmetry-breaking 
parameters, such as e. Such schemes are able to account 
for both a hierarchical pattern of eigenvalues, and order 
unity, sometimes maximal, mixing. Mass matrices have 
also been proposed where precise numerical ratios of 
different entries lead to the desired hierarchy and mixing. 

In this letter we propose an alternative view. This 
new view selects the large angle MSW solution of the 
solar neutrino problem, which is preferred by the day 
to night time flux ratio at the 2a level ||. While the 
masses and mixings of the charged fermions certainly 
imply regulated, hierarchical mass matrices, we find 
the necessity for an ordered structure in the neutrino 



sector to be less obvious. Large mixing angles would 
result from a random, structureless matrix, and such 
large angles could be responsible for solar as well as 
atmospheric oscillations. Furthermore, in this case the 
hierarchy of Am 2 need only be an order of magnitude, 
much less extreme than for the charged fermions. We 
therefore propose that the underlying theory of nature 
has dynamics which produces a neutrino mass matrix 
which, from the viewpoint of the low energy effective 
theory, displays anarchy: all entries are comparable, no 
pattern or structure is easily discernable, and there are no 
special precise ratios between any entries. Certainly the 
form of this mass matrix is not governed by approximate 
flavor symmetries. 

There are four simple arguments against such a pro- 
posal 

• The neutrino sector exhibits a hierarchy with 
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solution, while Am 2 tm 
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• Reactor studies of V e at the CHOOZ experiment 
have indicated that mixing of v e in the 10 _3 eV 2 
channel is small 0, requiring at least one small 
angle, 

• Even though large mixing would typically be ex- 
pected from anarchy, maximal or near maximal 
mixing, as preferred by SuperKamiokande data, 
would be unlikely, 

• Ve^fj, and v T fall into doublets with ej,, [It, and T£, 
respectively, whose masses are extremely hierarchi- 
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By studying a sample of randomly generated neutrino 
mass matrices, we demonstrate that each of these argu- 
ments is weak, and that, even when taken together, the 
possibility of neutrino mass anarchy still appears quite 
plausible. 

2 We have performed an analysis of a sample of ran- 
dom neutrino matrices. We investigated three types of 
neutrino mass matrices: Majorana, Dirac and seesaw. 
For the Majorana type, we considered 3x3 symmetric 
matrices with 6 uncorrelated parameters. For the Dirac 
type, we considered 3x3 matrices with 9 uncorrelated 
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parameters. Lastly, for the seesaw type, we considered 
matrices of the form MdM^Mq where Mrr is of 
the former type and Mr> is of the latter. We ran one 
million sample matrices with independently generated 
elements, each with a uniform distribution in the interval 
[—1,1] for each matrix type: Dirac, Majorana and 
seesaw. 

To check the robustness of the analysis, we ran smaller 
sets using a distribution with the logarithm base ten 
uniformly distributed in the interval [—1/2,1/2] and 
with random sign. We further checked both of these 
distributions but with a phase uniformly distributed in 
[0, 2ir]. Introducing a logarithmic distribution and phases 
did not significantly affect our results (within a factor of 
two), and hence we discuss only matrices with a linear 
distribution and real entries. 

We make no claim that our distribution is somehow 
physical, nor do we make strong quantitative claims 
about the confidence intervals of various parameters. 
However, if the basic prejudices against anarchy fail in 
these simple distributions, we see no reason to cling to 
them. 

In each case we generated a random neutrino mass 
matrix, which we diagonalized with a matrix U. We then 
investigated the following quantities: 

R = Am 2 2 /Am 2 3 , (1) 
Sc =4|t/ e3 | 2 (l-|(7e3| 2 ), (2) 

Satm=MU^\ 2 (l-\U^\ 2 ), (3) 

s Q =4\U e2 \ 2 \U el \ 2 , (4) 

where Amf 2 is the smallest splitting and Am| 3 is the 
next largest splitting. What ranges of values for these 
parameters should we demand from our matrices? We 
could require they lie within the experimentally preferred 
region. However, as experiments improve and these 
regions contract, the probability that a random matrix 
will satisfy this goes to zero. Thus we are instead 
interested in mass matrices that satisfy certain qualitative 
properties. For our numerical study we select these 
properties by the specific cuts 

• R < 1/10 to achieve a large hierarchy in the Am 2 . 

• sc < 0.15 to enforce small v e mixing through this 
Am 2 . 

• Satm > 0.5 for large atmospheric mixing. 

• Sq > 0.5 for large solar mixing. 

The results of subjecting our 10 6 sample matrices, 
of Dirac, Majorana and seesaw types, to all possible 
combinations of these cuts is shown in Table First 
consider making a single cut. As expected, for all types 
of matrices, a large percentage (from 18% to 21%) of the 
random matrices pass the large mixing angle solar cut, 



Dirac 


no cuts 


Satm, 


s© 


Satm + Sq 


no cuts 


1,000,000 


671,701 


184,128 


135,782 


sc 


145,000 


97,027 


66,311 


45,810 


R 


106,771 


78,303 


17,538 


14,269 


s c +R 


12,077 


9,067 


5,656 


4,375 



Majorana 


no cuts 


Satm 


Sq 


Satm + Sq 


no cuts 


1,000,000 


709,076 


200,987 


164,198 


sc 


121,129 


91,269 


70,350 


56,391 


R 


200,452 


149,140 


37,238 


31,708 


sc+R 


21,414 


16,507 


12,133 


10,027 



seesaw 


no cuts 


Satm 


sq 


Satm + Sq 


no cuts 


1,000,000 


594,823 


210,727 


133,800 


sc 


186,684 


101,665 


86,511 


49,787 


R 


643,394 


390,043 


132,649 


86,302 


sc+R 


115,614 


64,558 


53,430 


31,547 



TABLE I. Mass matrices satisfying various sets of cuts for 
the real linear Dirac, Majorana and seesaw scenarios. 
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FIG. 1. The distribution of Am /Am;; tm for Dirac (solid) 
Majorana (dot-dashed) and seesaw (dashed) scenarios. 

and similarly for the large mixing angle atmospheric cut 
(from 59% to 71%). Much more surprising, and contrary 
to conventional wisdom, is the relatively large percentage 
passing the individual cuts for R (from 10% to 64%) and 
for s c (from 12% to 18%). The distribution for R is 
shown in Figure |l|. Naively, one might expect that this 
would peak at R — 1, which is largely the case for Dirac 
matrices, although with a wide peak. In the Majorana 
case there is an appreciable fraction (~ 20%) that have 
a splitting R < 1/10, while in the seesaw scenario the 
majority of cases (~ 64%) have a splitting R < 1/10 — 
it is not at all unusual to generate a large hierarchy. 

We can understand this simply: first a splitting of a 
factor of 10 in the Am 2 's corresponds to only a factor 
of 3 in the masses themselves if they happen to be 
hierarchically arranged. Secondly, in the seesaw scenario, 
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FIG. 2. Plots of the normalized, binned distributions of 
Satm for Dirac, Majorana and seesaw cases. Contrary to 
intuition, the distributions actually peak at large s a tm- 

taking the product of three matrices spreads the Am 2 
distribution over a wide range. 

While one would expect random matrices to typically 
give large atmospheric mixing, is it plausible that they 
would give near-maximal mixing, as required by the Su- 
perKamiokande data? In Figure |^ we show distributions 
of Satm, which actually peak in the 0.95 < s a tm < 1-0 bin. 
We conclude that it is not necessary to impose a precise 
order on the mass matrix to achieve this near-maximal 
mixing. Finally, we consider correlations between the 
various cuts. For example, could it be that the cuts on R 
and sc selectively pass matrices which accidentally have 
a hierarchical structure, such that s a tm and Sq are also 
small in these cases? From Table | we see that there 
is little correlation of s atm with sc or R: the fraction 
of matrices passing the s a tm cut is relatively insensitive 
to whether or not the sc or R cuts have been applied. 
However, there is an important anticorrelation between 
Sq and sc cuts; for example, in the seesaw case roughly 
half of the matrices satisfying the sc cut satisfy the 
Sq cut, compared with 20% of the original set. This 
anticorrelation is shown in more detail in Figure [|, which 
illustrates how the sc cut serves to produce a peak at 
large mixing angle in the s Q distribution. 

For random matrices we expect the quantity 

s c + s Q = A(\U el U e2 \ 2 + \U el U e3 \ 2 + \U e2 U e3 \ 2 ) (5) 

to be large, since otherwise v e would have to be closely 
aligned with one of the mass eigenstates. Hence, when 
we select matrices where sc happens to be small, we are 
selecting ones where s© is expected to be large. 
3 We have argued that the neutrino mass matrix may fol- 
low from complete anarchy, however the electron, muon, 
tau mass hierarchies imply that the charged fermion mass 
matrix has considerable order and regularity. What is 
the origin for this difference? The only answer which we 
find plausible is that the lepton doublets, (fj, appear 



cuts 


none 


Satm 


s© 


Satm + Sq 


none 


1,000,000 


537,936 


221,785 


126,914 


sc 


222,389 


102,178 


99,050 


50,277 


R 


643,127 


345,427 


142,789 


81,511 


s c + R 


143,713 


65,875 


63,988 


32,435 



TABLE II. Mass matrices satisfying various sets of cuts for 
the real linear seesaw scenario, with additional mixing from 
the charged lepton sector. 



randomly in mass operators, while the lepton singlets, In, 
appear in an orderly way, for example, regulated by an 
approximate flavor symmetry. This idea is particularly 
attractive in SU(5) grand unified theories where only the 
10-plets of matter feel the approximate flavor symmetry, 
explaining why the mass hierarchy in the up quark sector 
is roughly the square of that in the down quark and 
charged lepton sectors. Hence we consider a charged 
lepton mass matrix of the form 

/A e \ 
Mi = Mi A M (6) 

\ X T J 

where A eiMiT are small flavor symmetry breaking param- 
eters of order the corresponding Yukawa couplings, while 
Mi is a matrix with randomly generated entries. We 
generated one million neutrino mass matrices and one 
million lepton mass matrices, and provide results for the 
mixing matrix U = Uju v , where U v and Ui are the 
unitary transformations on v\ and i; which diagonalize 
the neutrino and charged lepton mass matrices. We find 
that the additional mixing from the charged leptons does 
not substantially alter any of our conclusions - this is 
illustrated for the case of seesaw matrices in Table O. 
The mixing of charged leptons obviously cannot affect R, 
but it is surprising that the distributions for the mixings 
s a tm,Q,c are not substantially changed. 
4 All neutrino mass matrices proposed for atmospheric 
and solar neutrino oscillations have a highly ordered 
form. In contrast, we have proposed that the mass 
matrix appears random, with all entries comparable in 
size and no precise relations between entries. We have 
shown, especially in the case of seesaw matrices, that not 
only are large mixing angles for solar and atmospheric 
oscillations expected, but Am Q « 0.lAm 2 tm , giving an 
excellent match to the large angle solar MSW oscillations, 
as preferred at the 2a level in the day/night flux ratio. 
In a sample of a million random seesaw matrices, 40% 
have such mass ratios and a large atmospheric mixing. 
Of these, about 10% also have large solar mixing while 
having small v e disappearance at reactor experiments. 
Random neutrino mass matrices produce a narrow peak 
in atmospheric oscillations around the observationally 
preferred case of maximal mixing. In contrast to flavor 



3 



6 -__ 
5 - 
4 - 
3 - 
2 - 



without Sq cut 
with s@ cut 



(a) Dirac 




without 8q cut 
with Sq cut 



(c) seesaw 



FIG. 3. Plots of the normalized, binned distributions of sq for Dirac (a), Majorana (b), and seesaw (c) cases. The distribution 
after imposing the sc cut (solid) shows a greater preference for large sq compared with the original distribution (dashed). 



symmetry models, there is no reason to expect U e s is 
particularly small, and long baseline experiments which 
probe Am 2 atm , such as K2K and MINOS, will likely see 
large signals in V e appearance. If Am^ tm is at the lower 
edge of the current Superkamiokandc limit, this could be 
seen at a future extreme long baseline experiment with a 
muon source. Furthermore, in this scheme Anig is large 
enough to be probed at KamLAND, which will measure 



large v e disappearance. 
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